Structural Properties and Relative Stability of (Meta) Stable Ordered, 
Partially-ordered and Disordered Al-Li Alloy Phases 
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We resolve issues that have plagued reliable prediction of relative phase stability for solid-solutions 
and compounds. Due to its commercially important phase diagram, we showcase Al-Li system 
because historically density-functional theory (DFT) results show large scatter and limited success 
in predicting the structural properties and stability of solid-solutions relative to ordered compounds. 
Using recent advances in an optimal basis-set representation of the topology of electronic charge 
density (and, hence, atomic size), we present DFT results that agree reasonably well with all known 
experimental data for the structural properties and formation energies of ordered, off-stoichiometric 
partially-ordered and disordered alloys, opening the way for reliable study in complex alloys. 

PACS numbers: 71.15.-m, 71.23.-k, 61.66.Dk, 64.60.Cn 
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I. INTRODUCTION 

Predicting alloy composition-temperature {x vs. T) 
phase diagrams has important practical implications, and 
much progress has been made in the last decade via ah- 
initio study of phase stability. Such predictions are a 
sensitive test as the relative stability of metallic alloys 
depend on small enthalpy differences. Notably, the order- 
disorder temperature of a superstructure is related to the 
difference between the formation enthalpy of the com- 
pound Al,„Li„ (x = 71/(71 + 77i)) and mixing enthalpj*^ 
of disordered Ali^^jLij; alloy at the same composition 
Historically, DFT results were sensitive to basis sets and 
approximations?^ especially comparing ordering and mix- 
ing enthalpies from different codesi^ In particular, repre- 
sentations of charge densities are critical, especially for 
solid-solutions, where former "bad" results were not due, 
as often proposed, to the approximation to the disor- 
dered statei^ How the disordered phases is addressed re- 
mains an ongoing issue, such as the longstanding dis- 
crepancies between approximations to partially-ordered 
(off-stoichiometric and thermal antisites^) and disordered 
phases, e.g., the coherent-potential approximatioi*^^— 
(CPA), special-quasi random structures^i^iiiS (SQS), and 
cluster expansions^iii"— (CE). 

Al-Li alloys, with their unusual elastic and structural 
properties and commercial importance, have been the fo- 
cus of many theoretical and experimental studies. With 
doping to initiate precipitate growth in grains or grain 
boundaries, they are suitable materials for aerospace and 
automobile applications, due to their low density, high 
elastic modulus, and high strength-to-wcight ratio. Be- 
tween the stable endpoints of Al-Al (fee) and A2-Li 
(bcc), numerous structural phases exist or competcr^ 
e.g., the stable (3 AlLi (B32), rhombohedral Al2Li3 
(C33), monoclinic Al4Li9 (B2/m) and metastable (off- 
stoichiometric) (J'-AlaLi (LI2). The /3 AlLi, for example, 
is a promising candidate as an anodic material in high- 
energy density batteries; J', which appears in the mis- 
cibility gap between Al and /3, is used to precipitation- 



harden commercial alloys J^"— Substituting Li for Al not 
only makes the alloy less dense but increases unexpect- 
edly the elastic modulii^ even though the Young's mod- 
ulus of Li (14 GPa) is seven times smaller than that of 
Al (91 GPa). Also, the valence density leads to a bulk 
modulus of Al (83 GPa) five times larger than that of 
Li (15 GPa), giving rise to a sensitivity to basis sets if 
density is not represented propcrlyi^ 

Comparing reliably properties of all competing phases 
in alloys is critical, and Al-Li is a sensitive, sufficiently 
complex, and yet unresolved case. While there are some 
successes in describing the relative stability of the Al-Li 
ordered phases^"— a reliable description for the disor- 
dered and partially-ordered phases is lacking. For Al- 
Li disordered phases, the virtual crystal approximation 
(VCA)^ CPA^ and CE methodsii^ have been used to 
describe the x-dependcnt equilibrium volume and forma- 
tion enthalpies, but quantitative prediction of the lattice 
constants and mixing enthalpy remains a problem. 

Using an optimal site-centered basis for density and 
potentials)^ we present ab initio DFT calculations that 
compare reliably the phase stability (formation en- 
thalpies) of competing ordered, partially-ordered, and 
disordered phases, and quantitatively reproduce (with- 
out adjustable parameters) the unusual alloying effects, 
including the lattice constants, i.e., a vs. x, ill-described 
in the past. We also estimate the impurity formation en- 
ergies at a; ~ and 1 (i.e., the solution enthalpies) of the 
Al-Ali 

— a;Lia;: & Sensitive quantity due to differing elec- 
tronic nature of Al and Li. For general configurations, we 
describe more properly the topology of electronic charge 
density, and hence atomic size, charge and alloying ef- 
fects, especially in random alloys ^''^^ crucial in Al-Li for 
a reliable investigation of (meta)stable phases<^ We also 
resolve longstanding disagreements between various esti- 
mates of disordered energetics, e.g., CPAii^ and SQSi^ To 
showcase the predictive accuracy, we compare and con- 
trast our results to experimental data and those from lin- 
ear augmented plane wave (L AP W) augmented spher- 
ical wave (ASW),^ linear muffin-tin orbital (LMTO) 
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methods, atomic sphere approximation (ASA). 



II. COMPUTATIONAL DETAILS 

We use a Korringa-Kohn-Rostoker (KKR) Coher- 
ent Potential Approximation (CPA) code^ based on a 
weighted Voronoi polyhcdra (VP)^ basis defined from 
saddle-point radii (SPR) in the charge densityi^ With this 
optimal basis-set, a proper representation of the topology 
of charge density and hence atomic sizes yields energy 
differences insensitive to basis-set L-truncation<^i^ 

We include s-, p-, d- and /-symmetries in the KKR 
Green's functions spherical harmonic basis, which is trun- 
cated at Lniax=3, where L= {l^m). Energy integra- 
tion of Green's functions use a complex- energy Gauss- 
Chcbyshcv semicircular contour with 18 points. The Bril- 
lioun zone integrations use Monkhorst and Pack^l special 
k-point method using a 20 x 20 x 20 mesh. We use the 
von-Barth-Hedin2S local density approximation (LDA) 
as parameterized by Moruzzi, Janak and Williamsi^ 
For random alloys, the screened-CPA^ is used to in- 
corporate metallic screening from charge correlations in 
the local chemical environment. More details are given 
elsewhere i^i^i 

The potential zero, i.e., muffin-tin zero Vq, can dramat- 
ically affect stability prediction for spherical potentials. 
We use a variational definition (X stands for VP or ASA), 



(1) 



with a sum over all sites in a unit cell. We spherically av- 
erage functions over the solid angle (dfi = d9d(j)) within 
r| < -Res J the maximal region required in a site-centered 
method. We then accurately integrate over the intersti- 
tial region of arbitrary VP (full shape) via isoparametric 
integration^^ which is fast with error controllable to ma- 
chine precision. 

The total energies can be evaluate using weighted 
VPs, denoted by KKR-CPA(VP), or weighted- VPs ap- 
proximated by unequal ASA spheres, denoted by KKR- 
CPA(ASA). For comparison, we also provide equal ASA 
sphere results - often used in other ASA-bascd codes. 
Definition ([T]) yields kinetic energies that approach those 
of LAPWfiS. For both ordered and disordered (i.e., CPA 
and SQS) results, we find significantly improved predic- 
tions using optimal SPR basis within each VP, in com- 
bination with the Vq^, even if the remainder of the cal- 
culations is based on ASA, as exemplified for Al-Li. The 
CPA is applicable to an arbitrary x (both solid-solutions 
and partially-ordered states^), and configurational aver- 
ages can be enlarged via the DCA<^ 
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FIG. 1. (Color online) (Bottom) AiJ/ (meV/atom) vs. x 
of Al Ah-iLia:, and various fee- and bee-based compounds. 
KKR (LAPW-iS) results are filled (hollow) symbols; B32 
experimenli^ (*); Al phase (solid curve); groundstate hull 
(dashed lines); and Al-rieh (off-stoiehiometric) LI2 (pluses). 
Impurity formation energy Ai5imp (meV/atom) from a limit 
of a fit to AE/ (middle) Al in Li (0.9 < a; < 1.0) and (top) 
Li in Al (0 < a: < 0.1), values indicated in figure. 



III. RESULTS 



We report formation and mixing enthalpies (Ai?/), im- 
purity formation (or solution) energies (AiS^^p), and lat- 
tice constant (a) changes versus composition [x). We 
address the critical charge representation in electronic- 
structure methods to make predictions from various the- 
oretical techniques more consistent. For Al-Li, compar- 
ing to experimental data should be done with care, as the 
literature contains large scatter in various values. It is 
also well known that there is a Li deficiency that extends 
0.3 mm into the bulk from the surface in Al-Li alloys<^ 
Hence, experiments affected by this deficient region will 
not provide reliable data versus x. 
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TABLE I. AEf (in meV/atom) for (dis) ordered alloys with 
equal and weighted-SPR-VP spheres relative to elements in 
same phase, using the variational vq, Eq. IT}. CPA and SQS 
results are from this work. CC denotes combined corrections. 
Other theoretical and experimental results are also given. 



System 


Method 


Equal sphere 


SPR sphere 














A1-A2 


KKR-ASA(VP) 






—35.5 




Al 


LAPWii 






-62.6 






EMTOS^ 






-39.4 






KKR-ASA 






-19.0 




LI2/A1 


KKR-CPA(VP) 


—78.3 


-4.60 


-122.1 


-85.8 


AlsLi 


KKR-SQS-16(VP) 




-33.5 




-64.5 




LAPWii 






-113.0 






JjjiVl L v.^^" 






-111.0 






ASW-ASA2a 






- 109.0" 






KKR-CPA(ASA) 


-47.7 


-53.9 


-132.8 


-91.1 




KKR-SQS-16(ASA) 




-55.1 




-55.0 




LMTO-ASA-CO^l 


—75.6 


n/a 


— 133.7'' 






LMTO-CPA-ASA2i 


—57.4 


+91.1 


— 282. 9^" 


— 264.7" 


Llo/Al 


KKR-CPA(VP) 


-71.8 


-54.2 


-142.1 


-115.7 


AlLi 


KKR-SQS-8(VP) 




-55.7 




-116.0 










— 139.5 












-156.8 






KKR-CPA(ASA) 


—58.3 


—32.3 


— 159.6 


— 127.4 




KKR-SQS-8(ASA) 




-47.2 




-98.4 




LMT0-CPA-LDA2^ 








-102.0 






— 187.8 


—36.3 


—264.0 


-143.0 


AlLi 


KKR-SQS-8(VP) 




-65.7 




-121.0 




Expt^ 






-251.9 






LAPWi^ 






—232.3 






EMT02^ 






—295.0 






KKR-CPA(ASA) 


-179.7 


-39.4 


-279.6 


-142.9 




KKR-SQS-8(ASA) 




—67.2 




-152.4 


B2 AlLi 


KKR-CPA(VP) 


— 118.2 




— 159.5 






LAPWlii 






-142.1 






EMT02^ 






-189.6 






KKR-CPA(ASA) 


-109.2 




-138.5 






LMTO-ASA-CGii 






-141.3 




LI2/AI 


KKR-CPA(VP) 






-70.3 


-92.7 


AlLia 


KKR-SQS-16(VP) 








-78.3 




LAPWii 






-68.0 






EMT02^ 






-91.3 




A1-A2 


KKR-ASA(VP) 






+8.2 




Li 


LAPWiiiiS 






-6.8 






EMT02^ 






-3.9 






KKR-ASA 






+4.1 





° Adjusted spheres via pseudopotential theory, see text. 
' Charge-neutral, adjusted spheres with CC and Lmai ~ 2. 
Charge-neutral, adjusted spheres with hmax = 2. 



A. (Dis)Ordered Formation Energies 

Our results will clearly demonstrate that KKR- 
CPA(VP) correctly predicts the stability of ordered, dis- 
ordered and partially-ordered phase all within the same 



fast codei^ In Figure [T] (bottom panel) we show the Ai?^ 
versus x for ordered AlmLi„ (x = n/{n + m)) and dis- 
ordered Ali_2;Li^ alloys. Selected values of the results 
arc provided in Table HI and compared with LAPW and 
available experimental data, as well as ASA results, such 
as from linear and exact^ MTO (EMTO) methods. 

Importantly, using an SPR optimal basis and Eq. ([1}, 
the KKR-CPA(VP) results for compounds are now in 
excellent agreement with LAPW resultsi^ and mea- 
sured datai^ For example, the measured IS.EJ^'^ is 
-252 ± 10 meV and our optimal-basis KKR(VP) re- 
suh is -264 meV, while LAPW finds -232 meV. For 
metastable LIq AlLi, Ai;^^" is -142 meV for KKR(VP) 
and -140 rneV for LAPW. AEf^ for metastable B2 
(CsCl) phase, predicted to exist at high pressurefi^ is 
-160 meV for KKR(VP) and -142 meV for LAPW. 

The KKR-CPA mixing enthalpy AHf^{x) for disor- 
dered fee Al-Li is shown by the solid (red) curve in 
Fig. [T] Generally, in phase diagram assessments, the 
solid-solution enthalpy of a fixed phase is fitted to a poly- 
nomial in X, such as, 

Air/i(a;) = ^L„„ x™(l-a;)" V (m + n) < 4 (2) 

where Lmn are the polynomial coefficients. For Al Al-Li, 
the asymmetry in the AH-j^{x) in Fig. [1] arises from the 
Madelung contributions to the total energy, which is pro- 
portional to (AQi)^, the effective charge on a site, which 
also has a corresponding asymmetrical form as a function 
of X. Therefore, site charge is an important quantity in 
Al-Li, and a proper description of charge density and the 
corresponding transfer effect is necessary, as will be dis- 
cussed below. 



1. (Off ) Stoichiometric Al-rich LI2 

The observed miscibilitjJ^ of Ll2-Al3Li with respect 
to Al-Al and B32-AlLi is reproduced. Indeed, the segre- 
gation line between pure Al and B32 sits at —132 meV at 
25% Li, whereas LI2 is —122 7TieV, i.e., barely unstable 
but a very low-energy excitation; and, concomitantly, the 
stability of Ll2-Al3Li relative to Al-Alo.75Lio.25 yields 
the observed LI2 precipitation within the a(Al) + /? two- 
phase fieldii^ LI2 remains metastable as free energies of 
the competing ordered phases change little because their 
order is near perfect whereas that of the Al-rich disor- 
dered phase decreases rapidly due to the increasing sol- 
ubility of Li and the increasing configurational entropy. 

This argument can be supported by looking at the en- 
ergetics for Al-rich, off-stoichiometric LI2, shown as plus 
symbol in the lower panel of Fig. [TJ where the excess Al 
goes only to a Li corner-site and the three Al face-sites 
remain fully occupied by Al. These data indicate that 
a partially-disordered arrangement on the Al-rich side is 
more stable than the Al phase, and, by comparing to the 
groundstate hull, slightly more favorable than Q!(A1) + fi 
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alone. Based only on K enthalpy differences, we di- 
rectly predict a higher stability of off-stoichiometric LI2 
between 8-20%Li with a peak at ■~15%, as observedi^ for 
the S'^Al^hi) precipitates, where it undergoes a transfor- 
mation to a(Al) -H(5(AlLi) state at '-618 KJI At the 15% 
peak, we estimated the maximum order-disorder temper- 
ature from Ai?^^ ^ ^^^^"^ as 581 K. Also, at perfect sto- 

ichiometry ASf^'^^" is 42 meV (or 488 K) and it com- 
pares well with the 522 K obtained via a much more 
involved CE fits using LAPW results,— both of which 
should be compared to the measured eutectic tempera- 
ture (869 K) showing the liquid phase is more stable than 
the homogeneous solid-solution. 



2. CPA vs. SQS 

For completeness, we compare the CPA to the SQS, 
within the same code and approximations. The SQS is 
an ordered cell meant to approximate the two-site atomic 
correlations of a random alloy to a small cutoff, but it 
does not include configurational averaging as in the CPA, 
except as offered by the limited number of inequivalent 
sites within the ordered local environments. In Table U 
we compared our SQS and CPA results at 25% and 75% 
(a 16-atom celU^) and at 50% (an average of two 8-atom 
configurations^, one state having large positive formation 
energy and the other negative, as is typical). They agree 
from 2 — 20% depending on composition x, but only if the 
SPR basis is used for both, as can be verified in Table H] 
both CPA and SQS are sensitive to basis-representation 
but agree reasonably well within the optimal SPR basis, 
especially if larger SQS cells are used to improve the ef- 
fective configurational average (see FCC 50% case, e.g., 
in Table U where VP results shows much improved agree- 
ment). 



where Lmn are the restricted fitting coefficients. The 
slopes at each endpoint are very different from one an- 
other and cannot be reliably extracted by a polynomial 
fit for all X from Eq. as typically done. 

Using Eq. © and our KKR-CPA(VP) resuhs, 
our Ai;Lijj"^i is -244 meV/atom and A£;A|p'" is 
—82 meV/atom . A thermodynamic fit to experimen- 
tal data using CALPHAD yields^ -222 mcV/atom for 
AEl^l''^^ and -154 mcV/atom for AE^^^'^^K Our 
present results now agree much better with this one set 
of estimated experimental data, where notably only the 
Al-rich side has Al measured data. In the Li-rich end, 
we have ignored relaxation effects around the impurity, 
which leads to a decrease of 19 meV for the Ll2-AlLi3 
compound^ suggesting a potentially small increase in the 
magnitude of our AE^^™ ^\ 

Midownik^ has shown, however, that phase diagram 
fitting is not always able to predict structural energy 
differences accurately; hence, there is significant varia- 
tion in the literature. For example, for Li in Al and Al 
in Li, other CALPHAD fits yield —131 meV/atom and 
-1-972 meV/atomf^ respectively, the later being unphys- 
ical; and, another CALPHAD fit yields —954 mcV/atom 
and —130 meV/atom,— respectively, giving a phase dia- 
gram not in very good agreement with the experiment 1^ 

Other theory results are in quantitative disagreement 
to experiment and the present results; for example, pre- 
vious SQS resuhs for ASj^i^j" Ai _^qq „^cV/atom, 
and cluster expansions yields —353 TneV/atom, while for- 
mer equal-sphere CPA results give —358 meV/atomi^ All 
previous values were consistent, but none were obtained 
from proper impurity limit. Many of these discrepancies 
between results of the various electronic-structure based 
methods (CPA, CE, SQS, etc.) have been explained us- 
ing Al-Ag alloys as case study^ where charge issues were 
not significant; whereas for Al-Li the charge representa- 
tion is crucial, as also noted for "big-atom/small-atom" 
alloys,^ 



B. Impurity Formation Energies 

We now address the solution enthalpy (or impurity for- 
mation energy) , a more sensitive quantity to the approxi- 
mations used, because they are related to the slope of the 
formation energy in the impurity limits. In the top and 
middle panel of Fig.[Tl we show AEf vs. x for Al phase in 
the two extreme limits (0 < x < 0.1 and 0.9 < a; < 1.0). 
The solution enthalpy must be determined from Eq. ^ 
using two restricted fits in the limits, i.e., for B in A 
(0 < a; < 0.1) and vice-versa (0.9 <x < 1.0), 



AE. 



AE; 



3:^0(1) 



^i^-^W x™(l-x)^ 



(3) 



B in A 
imp 



AE; 



A in B 



x=0 



C. Lattice Constants 

The minimized a vs. x for Al Ali-^Li^^ (solid curve) 
are shown in Fig. [21 as well as for bulk LI2 AlaLi (square), 
with dilute limit slope added (dashed line) reproducing 
the well-known strong deviation to Vegard's rule (dot- 
dashed line). Although the atomic volume of Li is 20% 
larger than that of Al, addition of Li to Al initially 
causes a contraction of a vs. x, which we also find 
in quantitative agreement with experiment j ^^'^'^'^^ Var- 
ious experimental a's report e d^^'^"" — are indicated by 
symbols, including Al-rich Al and (^'(AlsLi) precipitates. 
First, our calculated a for pure Al at K is in agree- 
ment with experiment 1^ Our results do not include ther- 
mal expansion effects, which are in room temperature 
measurements, explaining the —0.7% difference in calcu- 
lated results. Experimentallyfi^ a for 5'{A\^1j\) at R.T. 
shrinks by 1% compared to pure Al, sec Fig. [21 and we 
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FIG. 2. (Color online) KKR-CPA(VP) a vs. x for Al- 
Ali-i^Lii; [solid (red) curve], compared to other data, includ- 
ing a^'pt at K.— a^pt from room temperature (RT) mea- 
surements [down^ and up^— triangles], aoxpt for ^'(AlsLi) 
precipitate at 523 K [*^]. Also shown are pseudo-potential 
results via first-order, perturbation theory [open circles].— 



TABLE II. Excess charges on the equal and weighted SPR 
of Li (first line) and Al (second line) atoms in (dis)ordered 
Al-Li. 



System 


Method 


Equal 


SPR 






^Qord 


AQ''" 




AQ''" 


LI2 AlsLi 


KKR-CPA(VP) 


-0.450 


-0.291 


-0.126 


-0.108 






+0.150 


+0.097 


+0.042 


+0.036 




KKR-CPA(ASA) 


-0.498 


-0.303 


-0.141 


-0.132 






-f0.166 


+0.101 


+0.047 


+0.044 


Llo AlLi 


KKR-CPA(VP) 


-0.101 


-0.069 


-0.032 


-0.026 






-fO.lOl 


+0.069 


+0.032 


+0.026 




KKR-CPA(ASA) 


—0.130 


—0.099 


—0.082 


—0.058 






+0.130 


+0.099 


+0.082 


+0.058 


B32 AlLi 


KKR-CPA(VP) 


-0.203 


-0.187 


-0.050 


-0.045 






+0.203 


+0.187 


+0.050 


+0.045 




KKR-CPA(ASA) 


-0.236 


-0.197 


-0.098 


-0.105 






+0.236 


+0.197 


+0.098 


+0.105 


B2 AlLi 


KKR-CPA(VP) 


-0.080 




-0.025 








+0.080 




+0.025 






KKR-CPA(ASA) 


-0.100 




-0.068 








+0.100 




+0.068 





find a reduction of 1.25% at K, in agreement with ex- 
periment (accounting for thermal expansion makes the 
agreement even closer). Finahy, a very sensitive mea- 
sure is the slope of the a vs. x curve near aluminium, 
which is measured^Sii^i^ to be -2.9 to -6.9 x 10"^ nm 
per at.%Li. The FLAPW impurity doping resultsii 
give —5.1 X 10^^ nm per at.%Li. We find it to be 
— 11 X 10^^ nm per at.%Li, within our CPA calculation. 



D. ASA versus VP 

For completeness, we provide the energetics from a 
wholly ASA implementation, including with Vq^^ in 
Eq. ^ evaluated using ASA spherical (not VP) inte- 
grals, denote as KKR-CPA(ASA). Moreover, to better 
understand the sensitivity of the effective charge and 
Madelung energy, the calculations are done in two differ- 
ent ways: (1) Standard equal atomic size (sphere or VP) 
for both Al and Li atoms (as done in many calculations), 
and (2) adjusting the sizes of each site by a weighted 
VP from the SPR found in a given environment 1^ From 
Table |T1 the SPR basis yields significantly improved 
energetics (especially compared to equal sphere ASA 
case) for both ordered and disordered phases, and elim- 
inates the basis-set dependence for energy differcncci^ 
For instance, the SPR KKR-ASA result for ordered 
Lls-AlgLi AS"^ is -133 meV/atom and it is in dra- 
matic contrast to the equal-sphere KKR-ASA result of 
—47.7 meV/atom. Improving the SPR- ASA basis using 
v^^ and VP integrations (i.e., KKR-CPA(VP)) yields 
AE^^^ = —122 meV/atom, in better agreement with 
-113 ?«eV/atom from LAPWJ^ 



E. Charge Representations 

The phase stability and the bonding characteristics of 
the Al-Li alloys can be understood by looking at distribu- 
tion of valence charge around each atom. A characteristic 
of this compound is that Li-atom redistributes some of 
its valence electron in between the Al bonds and the re- 
sultant strengthened Al bonds stabilize the compounds. 
The charge density topology in the Li-rich compounds is 
more crucial in stabilizing a particular phase. 

In Table ini the calculated local excess charges are pro- 
vided within ASA and weighted- VP, for some ordered 
and disordered phases. From the data, there is clearly 
a significant "charge transfer" from Li to Al for both 
ordered and disordered phases, as expected from elec- 
tronegativity where Al (1.61) is larger than Li (0.98). 
The excess (deficit) SPR- VP site charges from the KKR- 
CPA(VP), as expected, reflect a better charge neutral- 
ity than the SPR-ASA site charges more approximate 
KKR-CPA(ASA). The KKR-CPA(VP) thus provides an 
improved estimate of the chemical potential and energet- 
ics, as shown, as well as charge effect due to the more 
proper charge density representation. In a VP-based cal- 
culations, Madelung contribution to the total energy, as 
compared to an ASA, is minimized. 

For the optimal SPR basis, the inscribed sphere re- 
flect more appropriately the extent of the charge density 
on an atomic site, and hence a more reliable estimate of 
site charges; also the case for the unequal SPR-ASA re- 
sults, but the energetics are sensitive to VP integrations. 
For example, in AlLi for the SPR- VP basis, the excess 
(deficit) charge on Al (Li) remains almost unchanged for 
A2 and B32 phases (with only 0.005 e~ exchanged to Li 
from Al); whereas for the SPR-ASA sphere the charge 
remains small but is double that of the VP result. 
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For the equal-sphere ASA resuhs, however, there is 
an unphysical deficit charge on the smaher atom (Li) 
from the improperly described tails of the charge density 
of larger atom (Al), which have been arbitrarily cut off 
at the smaller radii (hence the dependence upon charge 
representation of the basis set we discussed) . The equal 
sphere results are unphysical, and the large charge ex- 
change significantly impacts the total energy results, see 
discussion in Section HVl which is also responsible for the 
enormous variation in past DFT results. 

IV. DISCUSSION 

Unlike other arbitrary (or unphysical) choices of charge 
density representation and basis set , the present SPR- 
based basis (weighted VP or ASA spheres) provides 
a unique, physics-based optimal representation of the 
charge topology for each atom type in a given environ- 
ment. This new basis reflects more properly the elec- 
tronegativity ("charge transfer"), reduces the overlap er- 
ror for ASA, and is valid for both ordered and disordered 
alloys, i.e., CPA and SQS approximations. For the sen- 
sitive case of Al-Li this has been demonstrated, which 
permits direct calculation of fully ordered, partially or- 
dered and disordered alloys, as recently shown for the 
quantum criticality in the doped, intermetallic NbFe2^ 
and magnetic-storage alloys <^ 

Previously, Masuda-Jindo and Terakurap2S applied 
pseudo-potential-based, first-order perturbation theory^ 
to reveal a contractions in a versus x for Al-rich solid solu- 
tions, similar to experiment, which deviates from simple 
Vegard's rule due to effects from bulk moduli and atomic 
volumes. They derived a relation for a mean ASA radii 
versus x for the alloy, which was then used to perform a 
KKR-ASW study2£ of solid-solution hardening and soft- 
ening using only ordered structures (Al7Li and AlsLi). 
For Ll2-Al3Li they found a decrease of a by 1.5% com- 
pared to pure Al. In light of our present results, their 
adjustments to the ASA radii better reflects the elec- 
tron density and "charge transfer" in agreement with our 
SPR-based VP and ASA. 

Recently, Korzhavyi et al^ performed LMTO-CPA- 
ASA calculations for the random alloy using charge cor- 
rections (similar to the scr-CPA) with an adjustable pa- 
rameter, and a posteriori correction for apparent charge 
transfer given by their choice of ASA spheres. Only then 
did they find that the a versus x and mixing enthalpy 
started to agree qualitatively with experiment. No such 
a posteriori choices and corrections are needed for our 
optimal basis, uniquely chosen before any calculationi^ 

Finally, in an attempt to correct the apparent "er- 
rors" of the CPA, Singh and Gonis^i proposed "charge- 
neutral" spheres, an arbitrary, unphysical and, as re- 
cently discussed)^ unnecessary assignment. Their for- 
mation energy changed from positive to negative when 
the spheres radii were adjusted to give charge-neutrality, 
with the latter giving a clearly incorrect formation en- 



thalpy of —283 (—265) meV for the ordered (disordered) 
AlaLi, as shown in Table HI For the ordered alloys, how- 
ever, they were able to add combined- corrections to the 
LMTO-ASA, which reduces ASA overlap errors due to 
the approximately 50% increase in the Al ASA radii. 
With combined-corrections the LI2 formation enthalpy 
improves to —134 meV; yet, they could not correct the 
CPA in a similar fashion. Such an approach, however, 
does not represent the CPA charge density correctlyi^ 
Simply put, for Al-Li, such a dramatic change in ASA 
radii to force charge neutrality arises solely because of the 
large difference in calculated bulk moduli (Al 0.72 Mbar 
vs. Li 0.12 Mbar) and volumes (Al 3.981 A and Li 
4.255 A). Hence, only for ASA, the Al sphere will enlarge 
at the expense of Li because the Li sphere is easily com- 
pressed (the ratio of Al to Li bulk moduli is 6), so the size 
of the Li sphere is expected to decrease much more than 
should be physically expected from the total density - see 
original motivations discussed in Ref. [13 and Ref.H. So, 
when the moduli are substantially different, this adjust- 
ment does not describe the proper density, and adjusting 
spheres to minimize energy, say, is unwarranted. 

Generally, however, using the saddle-point radii from 
the electronic density as the inscribed sphere of the VP 
(and using radical plane construction weighting to define 
full VP) offers a unique, physical and optimized repre- 
sentation, working for all ordered and disordered config- 
urational averaging approximations, including the CPA. 



V. SUMMARY 

For the commercially important Al-Li system, using 
a unique, physical and optimized representation of a 
site-centered charge density and potential in any local 
configuration, we accurately predicted, via KKR-CPA, 
the relative stability of all ordered, off-stoichiometric 
partially-ordered and disordered phases, including ther- 
mal antisitesi^ We also resolved long-standing discrepan- 
cies throughout the literature for the sensitivity of Al-Li 
energy estimates for disorder alloys (CPA, SQS, and CE), 
permitting accurate prediction of formation enthalpies, 
solution enthalpies and structural properties for general 
alloy configurations, all within a single code as desired 
for a reliable study of phase stability in complex alloys. 
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